Dusty waves on a starry sea: the mid-infrared view of M31 
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O ' ABSTRACT 
O 
(N 

Mid- infrared observations of the Andromeda galaxy, M31, obtained with the 
Q ! Infrared Array Camera on board the Spitzer Space Telescope are presented. The 

image mosaics cover areas of approximately 3? 7 x 1?6 and include the satellite 
galaxies M32 and NGC 205. The appearance of M31 varies dramatically in 
\ the different mid-infrared bands, from the smooth bulge and disk of the old 

stellar population seen at 3.6 /im to the well-known '10 kpc ring' dominating 
Q\ ', the 8 /im image. The similarity of the 3.6 fim and optical isophotes and nearly 

constant optical-mid-infrared color over the inner 400" confirms that there is 
no significant extinction at optical wavelengths in M31's bulge. The nuclear 
q \ colors indicate the presence of dust but not an infrared-bright active galactic 

"^h nucleus. The integrated 8 /im non-stellar luminosity implies a star formation 

rate of O.4M yr _1 , consistent with other indicators that show M31 to be a 
quiescent galaxy. 
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Introduction 



The proximity of the Andromeda galaxy, M31, has long made it a prime target for studies 
of galaxy structure and stellar populations. One tec hnique has been to isol ate particular 
populations: for example, studies of luminous stars ([Humphreys et al.lll990l ) showed that 
M31 has a lower massive star formation rate than M33 and the LMC. Global studies of 
particular co mponents have also b een valuable: H I maps provide estimates of the mass 
distrib ution (IBraun &: Thilkerll2004r) and have recently rey ealed a population of high- ve locity 



clouds (Thilker eta 



IRAS (IHabing et al 



2004). ISOPHOT flHaas et alJll998f ). MSX (IKraemer et al.ll2002h . and 



19841 ) observations showed a bright ring of infrared emission with radius 



10 kpc coincident with many of the H II regions. 

Mapping the entire disk of M31 at mid-infrared wavelengths allows bo th local and globa l 
studies of the galaxy. Observations with the Infrared Array Camera (IRAC; iFazio et al.ll2004l ) 
on the Spitzer Space Telescope simultaneously trace the dust in the spiral arms at 8 /im and 
the oldest stars in the disk and bulge at 3.6 and 4.5 /im without the complicating extinction 
or distance effects that make such studies in the Milky Way difficult. IRAC observations of 
M31 are complemented by deep data now available at many other wavelengths. They also 
complement Spitzer studies of other nearby galaxies. 

This paper presents an initial look at the IRAC observations of M31, focusing on the sur- 
face brightness profiles and extended emission. C ompanion papers discuss longer-wavelength 
MIPS observations of M31 dGordon et alj|2006h . IRAC and MIPS observations of the M31 
satellite galaxy NGC 205 (IMarleau et al.l 120061 ), and the implicat ions of the galaxy's mor- 
phology as seen in 8 fim non-stell ar emission ( Block et al.ll2006l ). A distance to M31 of 
783 kpc Jstanek fc GarnavichlE998l ) is assumed throughou t. All magnitudes are on the Vega 
system, using the calibration given by iReach et al.l (120051 ). 



2. Observations and Data Reduction 

The IRAC observations of M31 were taken as part of Spitzer General Observer program 
3126 in 2005 January and August^] Fifteen Astronomical Observation Requests (AORs) were 
used to map a region approximately 3?7x 1?6 (chosen to match the Spitzer/MIPS observations 
made as part of program ID 99), with an extension to the NW to include NGC 205. The 
central 1?6 x 0?4 was covered by three AORs, each having two 12-second frames per position. 



*A large solar proton event made most of the January data unusable; re-observations were made in 
August. 
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The outer regions were covered by two AORs each with two dithered 30-second frames per 
position. This mapping strategy ensured that observations of each position in the galaxy 
were separated by at least 2.5 hours, allowing efficient asteroid rejection in data processing. 
The complete dataset consists of 3000 individual images in each of the IRAC channels. 

Data reduction began with the Basic Calibrated Data (BCD) produced by versions 11 
(for the January data) or 12 (the August data) of the Spitzer Science Center (SSC) Pipeline. 
A 'delta dark' offset correction was applied to the 12- second 5.8 fim frames to correct for the 
first-frame effect, followed by use of the 'artifact corrector' software developed by S. Carey, 
which attempts to remove the electronic effects caused by bright stars. The remaining 
processing steps were carried out using the "IRAC_proc'fl software with the 2005 September 
5 version of the SSC MOPEX package. The 'overlap correction' necessary to make the sky 
background consistent between frames was computed and applied. Median-combined mosaics 
of the full dataset were created, then projected into the reference frame of each input BCD 
image to create 'source maps'. These source maps were used as the input 'uncertainty images' 
for creating the final photometric mosaics — this step prevents the outlier rejection used by 
MOPEX from rejecting pixels in bright sources and thereby biasing the photometry. The 
final mosaics were created with a pixel scale of Of! 86, sub-sampling the native IRAC pixel 
scale by a factor of y/2 to better sample the point spread function (FWHM ~ l'/9). The 
sky background was removed using the IRAF task IMSURFIT to model the background as 
a polynomial surface sampled by 'blank sky' regions near the edge of the mosaics. After 
subtraction, the mean image counts in small areas within these regions were computed, and 
the standard deviation of these means used to estimate the photometric uncertainties due 
to the background. 

The final step in mosaic production was photometric calibration. Althoug h calibration 



of po int source photometry varies slightly over the IRAC arrays' field of view (jReach et al. 



20051 ). we did not correct for this since doing so would have corrupted the extended emission. 
The calibration of IRAC is known to differ between point and extended sources, particularly 
at 5.8 and 8 fim where a significant fraction of point-source light is scattered on scales 
comparable to the array size. We used the results of recent work by T. Jarrett (priv. comm.) 
to correct from point to extended source calibration in the M31 mosaics by multiplying the 
images by (0.91, 0.94, 0.66 , 0.74) in the four IRAC bands, respectively [these values differ 



slightly from those given in lReach et al.l (120051 )] . The standard deviations of the background 



levels are about 0.04 MJy sr 1 in all bands except 5.8 /zm, where the background subtraction 
was more difficult, and the noise was 0.1 MJy sr _1 (lo"). These correspond to limiting surface 



2 IRAC_proc was developed by M. Schuster, M. Marengo and B. Patten at the Smithsonian Astrophysical 
Observatory. 
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brightness levels of 21.2, 20.9, 19.0, and 19.6 mag arcsec in the four IRAC bands. 

The final of mosaics are shown in Figure [TJ As with other recent works on IRAC 
observations of nearby galaxies we assumed that the 3.6 /im emission is purely due to stars, 
and used this to construct maps of the 'non-stellar' emission at 5.8 and 8 /im. The stellar 
emission at the two wavelen gths is assumed to correspond to 0.399 and 0.232 x the 3.6 (Ltm 
emission (IHelou et al.l 12004] ) . The 8 /im non-stellar image is shown in Figure [H 



3. Morphology and surface brightness profiles 



M31's appearance chang es dramatically over the IRAC bands, in a similar way t o other 
spiral galaxies such as M81 Jwillner et al.l 120041 ) and N GC 300 dlielou et al l bo04h. The 
3.6 /im view is similar to the i^-band image presented in Jarrett et al.l (120031 ): dominated 
by the bulge, with the spiral arms and disk relatively faint § At 8 //m, M31's bulge is much 
less prominent, and the spiral arms and the 10 kpc ring are the dom i nant features. The o uter 
(14 kpc) ring detected at far-infr ared wavelengths (|Haas et al.lll998t iGordon et al.ll2006l ) and 
tentatively detected with MSX (IKraemer et al.ll2002l ) is also seen in the 8 /im image. The 
8 /im image — par ticularly the non-ste llar version — looks extremely similar to the 24 /im 
image presented bvlGordon et al. ( 2006 ). showing the same split on the southwest side of the 
10 kpc ring that IGordon et all (booj T attribute to interactions with M32. The 'hot spots' 
seen near the nucleus at 160 /im are also prominent in the 8 /im non-stellar image, where 
they coincide with the ends of spiral arms. 

Disk images 3?3 x 0?9 were extracted from the mosaics and used to derive IRAC surface 
brightness profiles of M31. Peaks of bright stars, a l'-radius region around M32, and image 
regions with no coverage were all masked from the profile measurements. Profiles were 
measured using two methods: making 20"-wide cuts along the disk major and minor axes 
using the iraf task pvector, and fitting elliptical isophotes to the images using the ellipse 
task. Because ellipse can be unstable at low surface brightnesses, the isophote shape 
parameters (ellipticity, position angle, center) were constrained to be constant beyond the 
edge of the visible disk at semi-major axes > 0?88. The isophote fit at this distance has 
an ellipticity (e = 1 — b/a) of 0.73 and position angle of 38°. Profiles in the three longer- 
wavelength bands were measured on the ellipses fit to the 3.6 /im image, so that colors could 



3 The 2MASS Large Galaxy Atlas image of M31 is suspected of having an overly-faint disk due to problems 
with background subtraction (T. Jarrett, pers. comm unication). We do not compare the IRAC quantitative 
measuremen ts to those given by iJarrett et all ((21303), instead estimating M31's K — [3.6] color by M81's 



value of 0.3 ( Willner et al.l l2004L 
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be measured at the same spatial locations. 



The relative prominence of the M31 disk and bulge in the four IRAC bands is reflected 
in the major and minor axis profiles, shown in Figure [2J Out to a semi-major axis of ~ 400" 
(approximately the edge of M31's bulge), all four bands have very similar profiles. Beyond 
400", the 5.8 /xm and 8.0 /xm-band profiles decline more slowly and show more variability due 
to the clumpy nature of the dust emission. The 10 kpc ring is visible in all of the major-axis 
profiles; the 8.0 /im profile in particular appears highly asymmetric because of the split in 
the ring. The minor axis profiles are smoother because they sample mainly the bulge. Both 
major and minor axis p r ofiles are quite similar in shape to the optical profiles presented by 
Walterbos fc Kennicuttl (119871 ). indicating that dust extinction does not have a major effect 
on the optical profiles. The nucleus of M31 is smaller than the IRAC PSF, and does not 
separate from the bulge emission in the radial profile. 

The color profiles of M31 shown in Figure [3] also depict the disk/bulge separation. 
The median isophote colors within 400" are [3.6] — [4.5] = —0.14, [3.6] — [5.8] = 0.03 and 
[3.6] — [8.0] = 0. 13. The [3.6] — [4 .5] color is consistent with theoretical expectations for M0 
giant stars (see lPahre et al.ll2004f ) ; the [3.6] — [5.8] and [3.6] — [8.0] colors are slightly redder, 
presumably because of dust emission. At larger distances, the disk dominates the light and 
the colors become redder, particularly at the l ocation of th e 10 kpc ring. Comparing the 
3.6 /im profile to the r-band optical profile of iKentl (119871 ) shows that the r — [3.6] color 
is nearly constant over the central 400", with a luminosity- weighted mean of r — [3.6] = 
3.41 mag. This corresponds to r — K « 3.1, reasonably consistent with the predictions of 
Marastonl (120051 ). The colors of the M31 nucleus (measured in a 2"4-radius aperture, with the 
IRAC point source calibrations and aperture corrections applied) are [3.6] — [4.5] = —0.15, 
[3.6] — [5.8] = 0.19 and [3.6] — [8.0] = 0.51. The red colors at the longer wavelengths suggest 
that the nucleus has more dust th an the bulge, but are not the very red colors characteristic 



of an AGN (e.g. Istern et alJliopj ). 



Table [T] gives integrated flux densities and magnitudes of M31 as measured in the 
largest elliptical isophote (a =1?8). The uncertainties in these values include the effects 
of background variation and a 10% uncertainty in the absolute extended source calibration. 
The IRAC flux densities are consistent with t he COBE/DIRBE me asurements at 3.5 and 
4.9 fim (245 and 128 Jy with 20% uncertainties: lodenwald et al.l ll998h and the MSX A-band 
(8.3 /im) value (159 ± 32 Jy; Gordon et ali 2006). Ne ar-to-mid-infrared colors are another 
way to check our photometry: lAaronson et al.l (J1980J) measured H = 1.42 in a 1554" di- 



ameter circular aperture centered on M31. In the same apert ure the 3.6 /xm mag nitude is 
0.84, giving H - [3.6] = 0.58. Using M81's K s - [3.6] » 0.3 Jwillner et al.l I2004T) and the 
mean H — K s = 0.27 measured for early-type (SO/Sa/Sb) spirals by IJarrett et al.l (120031 ). 
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one would predict H — [3.6] = 0.57, in very good agreement with our measurement. 

A model light distribution consisting of an r 1 / 4 bulge and exponential disk were fit to 
the 3.6 jum semi-major axis surface brightness (SB) profile. The structural parameters fitted 
to the disk are sensitive to inclusion of the ring at semi-major axis lengths 2000 < a < 
3500 arcsec, while the bulge effective radius is sensitive to the inner fitti ng cutoff due to 
deviations from the idealized r 1 / 4 model at a < 100 arcsec (see iKentl 119831 ) . The sky value 
for the image was determined at isophotes roughly corresponding to 3.5 disk scale lengths. 
As there is still some galaxy emission at this radius, our sky value is likely over-estimated 
and is probably the dominant cause of the significant reduction in the SB at a > 4000 arcsec 
(smaller, but still significant, errors will be present in the profile at a < 4000 arcsec). We 
chose the fitting regions of 100 < a < 2000 and 3500 < a < 4000 arcsec to minimize these 
effects, and for consistency with previous work. 

The resulting model has (bulge, disk) half-light semi-major axis lengths of (10.3±0.4, 44. 7± 
0.6) arcmin, average ellipticities near to those scale lengths of (0.48 ±0.02, 0.68 ±0.04), total 
magnitudes (0.84 ± 0.13,0.58 ± 0.06), and RMS of 0.02 mag. The corresponding circular- 
ized effective radii are (7.5 ± 0.3,25.2 ± 0.7) arcmin. This disk major axis scale length of 



6 . 8±0.09 kpc compares we l 



with the longest wavelength (i?-band) value of 5.9 kpc measured 



by IWalterbos fc Kennicuttl (119881 ). The ratio of the bulge effective radius of 1.70 ± 0.07 kp c 



to disk scale length place M31 on the upper envelope of local spiral galaxies (Ide Jonglll996l ). 
The fitted bulge and disk models have a flux ratio B/D = 0.78 ± 0.11. This value is biased 
high relative to the true ratio, however, because the bulge model is t oo bright at a < 100 arc- 



sec an d the disk model is too faint at the location of the 10 kpc ring. IWalterbos fc Kennicutt 
(119881 ) calculate a similar value in the i?-band. 



The model-fitting allows an estimate of M31's total magnitude at 3.6 /zm through extrap- 
olation of the disk model to infinite radius. The magnitude within a < 4000 mag is 0.19 mag, 
the extrapolated disk from 2.5 scale lengths to infinity adds an additional 0.38 mag, and the 
sky subtraction errors on the SB profile within this isophote result in an additional 0.15 mag. 
The implied total magnitude of M31 is therefore —0.34 mag. We have chosen not to quote 
this model-dependent value in Table 1 for consistency with previous work and our measure- 
ments at other wavelengths, but it is consistent with the measurement in the 1?8 isophote 
given there when similar corrections are applied. The errors on the total magnitude of M31 
are dominated by the systematic error of sky subtraction on this large galaxy, and less so on 
the extrapolation of the disk emission. 

The 3.6 jum luminosity of M31 can be used to estimate t he ga laxy's stellar mass. Using 



a typical value of B — R = 1.5 fromlWa 



terbos fc Kennicuttl (Il987l ) and the color-dependent 



mass-to-light ratios of iBell fc de Jong! (120011 ). an estimated M/Lk for M31 is about 0.8 in 
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solar units. Using K — [3.6] = 0.3 and multiplying by M/Lk, the resulting stellar mass is 



Widrow et al. 



1.1 x 10 n M^. Recent dyn amical mass estimates of the M31 disk+bulge (e.g. 
20031 : iGeehan et aDbood ) are ~ 1O 11 M , in reasonable agreement with this derived stellar 
mass. Asymptotic giant branch ( AGB) stars can b e important contributors to the total near- 
infrared luminosi ties of galaxie s (lMarastonll2005l ). However, the most luminous AGB stars 
are predicted by iGroenewegenl (120061 ) to have much redder colors (0.5 < [3.6] — [4.5] < 1) 
than observed for M31, so we conclude that they do not dominate M31's 3.6 /im luminosity. 

The 8 /iin non-stellar emission is dominated by the 7.7 /im polycyclic aromatic hydro- 
carbon (PAH) band, which is a useful though imperfect star formation rate (SFR) indicator. 
It is therefore instructive to compare the M31 8 /im luminosity with other star formation 
indicators, suc h as Ha, r adio, and total infrared luminosities and with these values for 
other galaxies. IWu et al.l (120051 ) derived a calibration for SFR as a function of 8 /im non- 
stellar luminosity z/Lj,[8/im(dust)] (hereafter L 8 ) using correlations between L 8 and L(Ha) 
and L(1.4GHz). The IRAC 8 non-stell ar flux de n sity m easured for M31 corresponds 
to a luminosity of log(L 8 /L ) = 8.8. The IWu et al.l (120051 ) calibration yields an SFR of 
O.4M yr -1 . 

It is possible to compare M31's 8 /im-derived SFR with SF Rs from other ind icators, but 
since all such indicators have calibration uncertainties [e.g., the IWu et al.l (120051 ) 8 /im/SFR 
calibration is based on a small number of galaxies with a limited range of properties], we 
instead compare the obser ved properties directly . To estimate the Ha luminosity of M31, we 
convert the value given by lDevereux et al.l (119941 ) to the 780 kpc distance, multiply by 0.8 to 
account for [Nil] contamination, and multiply by 3.4 (as done by lWalterbos fc Braunlll994l ) 
to correct for extinction. The resulting L(Ha) = 2. 6 x 10 7 L^, which predicts log(L 8 /L ) = 
9.1. The 1.4 GHz radio flux density measured bv feeck et all Jl99sh. 3 34 x 10 20 W Hz" 1 , 
yields a predicted log(L 8 /L ) = 8.4. Figure 12 of iDale et al.l (120051 ) gives the ratio of 
vf v [8 / im ( dust ) ] / / ( 3 — 110 0/im) as a function of / I/ (70/xm)// I/ (160/im) for SINGS galaxies. 
From iGordon et all (bood ). M31 has /(70)//(160) = 0.12, with the corresponding /(3 - 
1100/im) = 2.57 x 10~ 10 W m~ 2 , yielding a predicted log(L 8 /L ) = 9.0. The Ha and far- 
infrared predictions are in reasonably good agreement, and are consistent with the observed 
8 /im luminosity given the scatter in the relationships. The 8 /im luminosity predicted from 
the radio flux is lower than observed, perhaps indicating that some extended emission was 
resolved out of the radio maps. All of the 'star-formation' luminosities are at the low end 
of the galaxy distribution. Despite its prominent star forming ring, M31 is a predominantly 
quiescent galaxy. 
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Fig. 1. — (Plate 1) Mid- infrared mosaics of M31 as seen in the IRAC bands. The grayscales 
are approximately logarithmic with intensity. From top to bottom: 3.6 /im, 4.5 /im, 5.8 fim, 
8 jttm, and 8 jum 'non-stellar' (produced by subtracting the expected stellar contribution to 
this band, a scaled version of the 3.6 /im image). All images are centered on the nucleus and 
cover an area of 3?0 x 0?87. NGC 205 was observed with IRAC but is not shown here (see 



Marleau et alJl200fl l. 
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Fig. 2. — Surface brightness profiles measured along the M31 disk major (left) and minor 
(right) axes. Profiles have been offset for clarity a s indicated in the legend. Th e mid- 
infrared profiles are compared to optical profiles from IWalterbos fc Kennicuttl (119871 ). which 
are plotted as flat in the center where no data were available. The upturn in the northwest 
minor axis profiles of IWalterbos fc Kennicuttl (119871 ) is due to NGC 205, not included in the 
IRAC images used here. 
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Fig. 3.— Mid-infrared surface brightness and color profiles of M31 derived from 
iraf/ellipse isophotes fitting. All data use the extended source calibration, as described 
in the text. '8.0ns' refers to the 'non-stellar' 8 /im image. (Top): IRAC surface brightness 
profiles. (Middle): IRAC color profile s from the sa me data. (Bottom): Optical- mid-in frared 
color profiles, using optical data from iKentl (119871 ) and IWalterbos fc Kennicuttl (119871 ). The 
large offset betw een colors using r and R is l ikely due to problems with the zero-point 
calibration of the IWalterbos fc Kennicuttl (119871 ) data. 
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Table 1. M31 mid- infrared flux densities and magnitudes 



A 


fu 




H 


[Jy] 




3.6 


259 ± 32 


0.09 ±0.13 


4.5 


144 ± 20 


0.24 ±0.15 


5.8 


190 ± 35 


-0.55 ±0.20 


8.0 


151 ±21 


-0.93 ±0.15 


8.0 (non-stellar) 


91 ±21 





Note. - - Measurements made in an ellip- 
tical aperture of semi- major axis a =1?8. 



